Anomalous critical behavior near the quantum critical point of a hole-doped La2Cu04 
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The Landau-Ginzburg- Wilson paradigm for critical phenomena is spectacularly successful when- 
ever the critical temperature is finite and all fluctuation modes, with characteristic energies much 
smaller than ksTc, obey classical statistics. In zero-temperature quantum critical phenomena, clas- 
sical thermal fluctuations are replaced by zero-point quantum fluctuations and quantum-mechanical 
generalization of the Landau-Ginzburg- Wilson paradigm has been a central topic in condensed- 
matter physics. In this neutron-scattering study on spin fluctuations near the quantum critical 
point induced by hole-doping in La2Cui-a:Lia:04 (0.04 < a; < 0.1), the phase boundary for quantum 
crossover expected from the generalized quantum theory for critical phenomena was observed for the 
first time. Furthermore, critical exponent and scaling function become anomalous near the quantum 
critical point, which has not been expected in current theories. 



Interesting phenomena, including high- Tc supercon- 
ductivity, occur in laminar cuprates below a tempera- 
ture T ^ J /ks ~ 1500 K, where J is the magnetic ex- 
change energy between the nearest-neighbor spins in the 
Cu02 plane. In this temperature range, quantum fluc- 
tuations become dominant and add one extra dimension, 
S^T = hc/kBT, where c is the spin-wave velocity, to the 
2-dimensional (2D) Cu02 square latticejll3| (Fig.[ll). At 
the zero temperature, a quantum phase transition from 
the antiferromagnetic (AF) phase to a paramagnetic one 
occurs at a critical hole doping, a.k.a. quantum criti- 
cal point (QCP), Xc, in the cuprates. Thus, for doping 
X > Xc on the paramagnetic side of the QCP, magnetic 
correlation length at T = 0, ^o, is finite. Generalizing the 
Landau-Ginzburg- Wilson paradigm for critical phenom- 
ena in classical systems to the quantum spin system in 
the (2+l)-dimcnsional space, the shorter of and sets 
the long-wavelength cutoff for spin fluctuations and is ex- 
pected to determine the universal magnetic phenomena 
with universal critical exponents^. The E/T scaling in 
spin dynamics, observed in diverse samples of hole-doped 
La2Cu04 using Ba, Sr or Li dopant, with superconduct- 
ing or insulating ground state, and with commensurate 
or incommensurate dynamic spin correlationsj3,IMSI3|' 
can be understood as the physical consequence of < Co 
at T > Tx = hc/kB^o- 

In this neutron scattering study on quantum spin dy- 
namics in hole-doped La2Cu04, we determine the dop- 
ing dependence of the crossover temperature Tx , below 
which the E/T scaling breaks down and crosses over to 
a constant-energy scaling, consistent with < Ct- The 
most interesting result of this study, however, is that the 
critical exponent a in the E/T scaling 3] changes from an 
expected a « 1 to an anomalous a ~ 0.65 when doping 
is reduced toward the AF QCP at Xc- The change of the 
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FIG. 1: (color) Phase diagrams for Sr and Li-doped La2Cu04 
where increasing hole concentration suppresses the antiferro- 
magnetic (AF) phase at the quantum critical point Xc and 
drives the Sr doped material superconducting (SC). The line 
filled area denotes spin-glass phase. We observe the crossover 
at Tx (red symbols) between the E/T scaling (QC) and the 
E/Tq scaling (QD) regimes in Li-doped La2Cu04. The inset 
illustrates the (2-|-l)-dimensional space formed by the Cu02 
plane (Cu; red circles, O: yellow circles) and the quantum 
^T- The blue hemisphere indicates the extent of the magnetic 
correlation length at zero temperature, ^o- 



critical exponent is inconsistent with the classical univer- 
sality of critical exponents. In analog to a recent example 
in heavy fermion metal "^l, the substantial reduction 
of the a exponent in the 2D spin system suggests extra 
physics in quantum critical phenomena which invalidates 
the classical Landau-Ginzburg- Wilson paradigm [loj. 

We choose Li-doped La2Cu04 single crystals for this 
study (Table I). Isotopically enriched ^Li (98.4%) was 
used to reduce neutron absorption of natural Li. The 
crystals have orthorhombic Cmca symmetry in the tern- 
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TABLE I: Doping x of La2Cui_iLia;04 and main experimen- 



X 


0.04 


0.06 


0.1 


Tx (K) 


35 


50 


30 


a 


0.65±0.06 


f»l 


«1 


b 


0.21±0.01 


0.18±0.01 


0.28±0.01 



perature ran ge o f this study and lattice parameters were 
reported in |T3|. Analogous to Sr-doped La2CuO4 0, 
Ha . the AF phase is suppressed rapidly by 3 per- 
cent of holes and a spin-glass transition occurs below 
10 K on both sides of the QCPjUEl (Fig-CJ. In ad- 
dition, similar spin dynamics following the E/T scaling 
have been observed ^, j6, 7\. To investigate the in- 
trinsic spin dynamics of the 2D quantum antifcrromag- 
net, the suppression of the d-wave superconductivit y by 
Li dopants as potential scatterers in the Cu02 plane [ig 
is crucial. This allows measurement of the momentum 
and energy-dependent dynamic magnetic structure fac- 
tor, S{q,E), at low temperature without interference 
by the superconducting gap as seen in Sr or Ba-doped 
La2Cu04!Ti,'l^. To achieve required energy resolution, 
the cold neutron triple-axis spectrometer IN14 at the In- 
stitut Laue-Langevin with Ef=5 meV (for a;=0.1), and 
SPINS at the National Institute of Standards and Tech- 
nology Center for Neutron Research with Ef=3.7 meV 
(for a;=0.04) were used. The scattered neutron beam was 
filtered using cooled Be or BeO to remove higher order 
neutrons. The dynamic spin correlations which can be 
resolved in this cold neutron inelastic scattering study 
do not participate in the spin- freezing process [3. We 
will report the minority spin-glass phase elsewhere. 

Dynamic spin correlations in Li-doped La2Cu04 are 
finite-sized fluctuating versions of the chessboard-like an- 
tiferromagnetic order in pure La2Cu04j23|. Intensity of 
^(q, E) concentrates in reciprocal space in rods inter- 
cepting the planar Brillouin zone of the Cu02 square 
lattice at commensurate Bragg points of the (tt, tt)- 
type. Scans at various E and T across the rod at 
Q=(7r,7r) are shown in the inset to Fig. |21 for Li dop- 
ing a;=0.1. Note that Q=(100) in the orthorhombic no- 
tation. The scans have been normalized by the peak 
intensities so that the invariance of the in-plane peak 
width in this work using coarse q resolution can be easily 
seen. Reflecting the two-dimensional nature of spin cor- 
relations, S{q,E) is flat in the out-of-plane direction^. 
Thus, / dqS{q, E) is simply proportional to the value of 
^(Q,^^) in the T and E range of this study, and the 
imaginary part of the local magnetic response function 
X"{E) = (1 - e--^/'^^^) / dc^S{q_,E) is proportional to 
(1 - e--^/'^^^) 5(Q, E). Measured S'(Q, E) as a func- 
tion of E and T is shown in Fig. |2t and b for x—Q.l 
and 0.04, respectively. Throughout this paper, a unique 
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FIG. 2: (color) Measured dynamic magnetic structure factor 
^(Q,^) as a function of E and T at Q=(100) in the Cmca 
orthorhombic notation which corresponds to the (vr, tt) point 
of the Cu02 square lattice, a, Li doping a;=0.1. b, 1=0.04. 
The blue curve at E=0 traces ^(Q, E) at the low energy limit 
as a function of temperature. Inset: normalized 5'(q, E) as a 
function of q at indicated E and T for a::=0.1. 



color is assigned to each temperature for data points in 
the figures. The peak at E=Q (dashed curves) contains 
incoherent, elastic and quasielastic scattering processes 
within the energy resolution, Ares, of the spectrometer 
and is not the subject of this work. For the resolvable 
dynamic signal, S'(q, E) with \E\ > A^es, background 
(straight green lines) is determined by measurement at q 
far away from the peak of >S'(q, E) (refer to the inset to 
Fig. 12), or a,t E ^ —ksT where the Bosc factor for spin 
fluctuation modes approaches zero. 

To demonstrate the E/T scaling in magnetic response 
function: 



X"{E)T- ^AfiE/bkBT), 



(1) 



where a is a critical exponent .3j , A a constant, /(y) a 
scaling function defined to peak at y=l, and b = Ep/ksT 
relates the peak energy of x"(i?) to thermal energy, ex- 
perimental data of x"{E) above Tx, renormalized by 
r~", are plotted semilogarithmically as a function of 
E/ksT for x=0.1 and 0.04 in Fig.|3t and b, respectively. 
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FIG. 3: (color) Scaling plots of the magnetic response func- 
tion, a, b, The E/T scaling at high temperature for x—0.1 
and 0.04, respectively. Inset in b shows the goodness-of-fit 
for the critical exponent a=0.65 for a::=0.04. c, d, Departure 
from the E/T scaling at low temperature. The solid line in 
a and c is Eq. Q using Eq. Q as the scaling function, in b 
and d using Eq. ^ as the scaling function. Color of symbols 
denotes temperature as in Fig. |5| 

For a::=0.1, identical to the case for x— a « 1 and 

the scaling function is 

hiy) = y/il + y'). (2) 

The peak position Ep/kBT=b = 0.28 ± 0.01 in Fig. ^ 
is larger than 5=0.18 for X—0.Q6. For Ba or Sr-doped 
La2Cu04, the exponent a has been determined only for 
x=O.U and it is also a = 0.94 ± 0.06 « 10. Thus for 
x>0.06, hole-doped La2Cu04 has a normal critical expo- 
nent a«l, consistent with a Curie's law for the staggered 
magnetic susceptibility. 

For Li doping a::=0.04 which is close to the AF QCP 
of a;c«0.03,ll], a fractional a=0.65 ± 0.06 is required for 
the data to collapse onto a single scaling curve (Fig. ISb) . 
Previously, in heavy fermion alloy CeCus.gAuo.i, which 
is also near an AF QCP, single-crystal neutron scat- 
tering data are collapsed onto the E/T scaling with 
a=O.74±O.10. The same procedure as in Ref. 8], which 
is unbiased with respect to the selection of f{y), is used 
to optimize a for a;=0.04, and the goodness-of-fit repre- 
sented by the chi-square statistical test is shown in the 
inset to Fig.|3j3. The spectrum in Fig.Os is broader than 
that in Fig. and cannot be described by Eq. (01. We 
adopt empirically the scaling function for x=0.04, 

f2iy)^sgn{y)y^/{l + yl), (3) 

where yi = \y\[{v - u)/^]^/" with m=0.77±0.02 and 
w=1.33±0.02 as the best fits. The parameter b = 0.21 ± 
O.Ol. Note that for \y\ <C 1, f2{y) ^ sgn{y)\y\^' is also 
anomalous and does not conform to the usual analytic 
form x"{E) ^ E at small E^. 



The heavy fermion metal CeCus.gAug.i {z=2 and di- 
mension d>2) and our a;=0.04 cuprate sample (z=l and 
d=2) belong to different universality classes with the ef- 
fective dimension (d + z) above and below the upper crit- 
ical dimension, respectively P, |^. However, they are 
both close to an AF QCP and share within error bars an 
anomalous a exponent substantially smaller than unity. 
In our study of doping dependence, furthermore, hole- 
doped La2Cu04 recovers the normal a « 1 and analytic 
x"{E) for X > 0.06. These changes in critical behavior 
are very unusual. Anomalous a exponent for the heavy 
fermion metal CeCus.gAuo.i is recently explained by in- 
cluding local degrees of freedom in the Kondo screen- 
ing, in addition to traditional long-wavelength magnetic 
fluctuations, in the antiferromagnetic quantum phase 
transition 9] . Closer to our case of doped cuprates near 
the AF QCP, currently, "deconfined" degrees of freedom 
associated with fractionalization of magnetic order pa- 
rameter in pure 2D antifcrromagnet are being explored 
theoretically 0. Like Ref. ^J, this theory goes beyond 
the Landau- Ginzburg- Wilson paradigm, which focuses 
on long-wavelength fluctuations of the order parameter in 
critical phenomena. Whether including these additional 
"deconfined" degrees of freedom in the critical theory for 
the QCP would, similar to the case for heavy fermion 
metals|2|, produce the anomalous a exponent and f2{y) 
reported here for cuprates is unknown. Nevertheless, our 
observation of anomalous critical phenomena near the 
QCP in another class of correlated electronic material 
suggests the generality of the phenomena, and would cer- 
tainly encourage the budding research interest on subtle 
quantum interference effects near a QCP. 

In Fig. O; and d, low temperature data are added to 
the scaling plots for a;=0.1 and 0.04, respectively. Since 
the peak position shifts to the right with decreasing tem- 
perature, the E/T scaling cannot be satisfied. Instead, 
as shown in Fig.Q^ and b, below Tx, the local magnetic 
response function x"{E) is nearly temperature indepen- 
dent. The curves in the figures are 

x"{E)^Bf{E/T^), (4) 

where i? is a constant, and the same pair of scaling 
functions, fiiy) and /2(y), are used in Eq. Q with 
ro=0.66±0.01 and 1.11±0.03 meV for a;=0.1 and 0.04, 
respectively. Note that there is no detectable gap in 
X' {E) in Fig. H for a;=0.1 and 0.04, as for x=0.06 re- 
ported previously Ji]. In some theoretical models of non- 
random quantum antiferromagnets, a quantum spin liq- 
uid has a gap, A ~ ?ic/^o, near E=0 in x"{E) below 
Tx |2l| . The gapless x"{E) observed in our samples is ob- 
tained in these theories with str ong damping of the spin 
liquid by doped holes, F « a |2^|23| . Partial filling of the 
Haldane gap in ID quantum spin liquid by doped holes 
has been reported lately In other theoretical mod- 
els of quantum antiferromagnets, several types of gapless 
spin liquids exist without impurity scatteringj25j. Clas- 
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FIG. 4: (color) Temperature invariance of the magnetic re- 
sponse function at T < Tx- a, The doping x—0.1. The sohd 
hne is Eq. Q using Eq. ^ as the scahng function, b, x=0.04:. 
The sohd hne is Eq. Q using Eq. @ as the scahng function. 
Color of symbols denotes temperature as in Fig. |21 



nomena to quantum critical phenomena in cuprates nat- 
urally explains the universal occurrence of the E/T scal- 
ing in hole-doped La2Cu04 and the crossover to the fi- 
nite magnetic-energy regime at low temperature '5,'^, the 
anomalous critical exponent and scaling function near the 
QCP reported here have yet to be understood, probably 
by going beyond the Landau- Ginzburg- Wilson paradigm 
and taking into account additional quantum interference 
effects existing near the zero temperature^, 
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sification of our observed 2D gapless spin liquid in hole- 
doped cuprates would be an interesting theoretical task. 

The crossover temperature, Tx, from the E/T scal- 
ing to the E/Tq scaling regime can be conveniently 
located by the peak in temperature dependence of 
5(Q, E) at the £'— »0 limitQ, since at high temperature, 
^(Q, £;^0) cx T-" from Eq. (HJ-® , and at low temper- 
ature 5(Q, E^O) cx T from Eq. ©-Q- The blue curves 
at E^O in Fig. |3i and b are S{Q,E^O) for a;==0.1 and 
0.04, respectively. The extracted Tx from Fig. |21 to- 
gether with previously determined value for a;=O.O60|, 
are shown in the phase diagram in Fig. ^ and are listed 
in Table I. Under the dome-shaped phase boundary, spin 
fluctuation energy, which would approach zero accord- 
ing to the E/T scaling, saturates at a finite Fq. Thus 
the crossover at Tx keeps spin fluctuations energetic at 
low temperature, which might be essential for the high 
value of Tc if superconductivity in cuprates is mediated 
magnetically by spin fluctuations. 

Experimental search for the E/T scaling was 
initially stimulated by the marginal Fermi liquid 
phenomenology 0, 0, Hg. The subsequent observation 
of the scaling in insulating sample supports the QCP 
as the common mechanism. The breakdown of the 
E/T scaling below a dome-shaped crossover boundary 
reported here excludes the existence of additional mag- 
netic QCPs in the investigated doping range. The 
gapped quantum spin liquid predicted below Tx for non- 
random 2D antiferromagnets does not survive in our 
doped cuprates. This observation disfavors those theo- 
ries which rely on the robustness of the gapped quantum 
spin liquid against doping for high- Tc superconductiv- 
ity in cuprates 27] . While the application of the classi- 
cal Landau-Ginzburg- Wilson paradigm for critical phe- 
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